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ABSTRACT 

We calculate for the first time the finite size corrections in the massive Thirring 
model. This is done by numerically solving the equations of periodic boundary 
conditions of the Bethe ansatz solution. It is found that the corresponding central 
charge extracted from the 1/L term is around 0.4 for the coupling constant of 
go — — f and decreases down to zero when g 0 = — This is quite different from 
the predicted central charge of the sine-Gordon model. 
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In two dimensional field theory, there is a remarkable correspondence between 
the fermionic and bosonic held theories. This was first recognized by Coleman 
[1], and he proved that the sine-Gordon held theory and the massive Thirring 
model are equivalent to each other in that the arbitrary order of the correlation 
functions turn out to be the same. 

Recently, however, Klassen and Melzer [2] argue that the equivalence between 
the sine-Gordon and the massive Thirring models may be violated at the finite 
size correction. They proved by using the perturbed conformal held theory that 
these two models are different in hnite-volume energy levels, for example. 

In this paper, we calculate the finite size corrections to the ground state 
energy. We solve numerically the equations of the periodic boundary condition 
in the Bethe ansatz solutions of the massive Thirring model [3-5]. The ground 
state energy can be expressed as 

7TC 

E v = E 0 L — — + ... (1) 

o L 

where L denotes the box size, c corresponds to a central charge at the massless 
limit [6,7]. 

The present calculation shows that the corresponding central charge c in the 
negative coupling constant regions (no bound states) is around 0.4 for g 0 = — | 
and that it becomes zero when go = —These values can be compared with 
those calculated for the sine-Gordon held theory [8,9]. The central charge for the 
sine-Gordon held theory with the massless limit can be expressed as 

i 6 

C = 1-t-— 

P{P + 1) 

where p is an integer and is related to the coupling constant go as 

9o =—(!--)■ (3) 

2 p 

In hg.l, we summarize the calculated central charge as the function of the coupling 
constant for the sine-Gordon model by Itoyama and Moxhay, and for the massive 
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Thirring model by the present calculations. One can see that the values of the 
central charge predicted for the two models are very different from each other. 

It is, however, not very clear to us whether this difference may be related to 
a possible violation of the equivalence between the sine-Gordon held theory and 
the massive Thirring model at the finite volume energy as suggested by Klassen 
and Melzer. 

Here, we briefly review the massive Thirring model whose lagrangian density 
can be written as [10] 

£ = ^(* 7 ^ - m ofin ( 4 ) 

with the fermion current j tl =: ^>7^ Choosing a basis where 75 is diagonal, 
we write the hamiltonian as 

c d (9 

H = dx +m 0 (^i^2 +V4'0i) + 20 O ^iV4^2^i • ( 5) 

J [ ox ox 

The hamiltonian eq.(5) can be diagonalized by the Bethe ansatz wave func¬ 
tions ^(xi, ...,Xn) with N particles 

T(a;i, ...,x N ) = exp(im 0 ^ Xi sinh Pi) II [1 + Pj)e{xi - Xj)] (6) 

where /% is related to the momentum k, and the energy E % of i-th particle as 


hi = mo sinh 

(7a) 

E t = mo cosh Pi. 

(76) 

where Pi’s are complex variables. 


e(x) is a step function and is defined as 


( — 1 x < 0 


e(x) = l 

(8) 

{ 1 x > 0. 
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( 9 ) 


A {fa, fa) is related to the phase shift function fafa — fa) as 


1 + i\(fa, fa) 
1 - i\(fa, fa) 


exp (i<f)(fa - fa)). 


The phase shift function 4>{fa ~ fa) can be explicitly written as 


<f>(fa - fa) = -2 tan -g 0 tanh -{fa - fa) 


( 10 ) 


From the definition of the rapidity variable fa's, one sees that for positive 
energy particles, fa's are real while for negative energy particles, fa takes the 
form in — on where ads are real. 

Since the Bethe ansatz wave functions diagonalize the hamiltonian, we de¬ 
mand that they satisfy the periodic boundary conditions (PBC) with the box 
length L [3], 

^(xi — 0) — ^(xi — L). ( 11 ) 


This leads to the following PBC equations, 


m 0 L sinh fa = 2nm - ^ 4>(fa - fa) (12) 

3 

where nds are integer. Here, we note that we cannot take the anti-periodic 
boundary condition since it does not reproduce the boson spectrum in the positive 
coupling constant regions [5]. 


The parameters we have here are the box length L and the particle number 
N. In this case, the density of the system p becomes 


P = 


N 

~L 


(13) 


Here, the system is fully characterized by the density p. 

We write the PBC equations for the vacuum which is filled with negative 
energy particles ( fa — in — a* ), 


sinh a* 


2nrii 

u 


2 

Lq 


y) tan 1 

j¥=i 


'I 

.2 


go tanh (aj o.j j , 


(14) 
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where n t = 0, ±1, ±2,±1V 0 with N 0 = |(7V — 1) and L 0 = rn 0 L. 
In this case, the vacuum energy E v can be written as 

N 0 

E v — — ^2 m o cosh a t . 

i=—No 


( 15 ) 


In this paper, we have carried out the numerical calculations of the PBC 
equations. The numerical method to solve the PBC equations is explained in 
detail in ref. [5]. 

Now, the calculated vacuum energy can be parametrized as 


E V = E 0 L-^ f^ + ... 

o L 


(16) 


where c(p 0 ) corresponds to the central charge at the massless limit. In what 
follows, we call this c(po) as the central charge even though we are solving the 
massive field theory. It should be noted that the first term in eq.(16) can be 
evaluated analytically by taking the thermodynamic limit [3]. 

Since we can vary the values of L and N, we obtain the corresponding central 
charge c(go). Although we have still rather small particle number (.N ~ 10000), 
we believe that the values extracted for the central charge must be reasonably 
reliable. 

Now, we want to obtain the central charge c(go) at the field theory limit 
p —*• oo. In fig. 2, we show the calculated central charge c(go) as the function of 
the effective density p 0 = It is quite interesting to observe that the calculated 
central charge can be well parametrized by the following simple formula [11], 

c(go) = A + B exp ( 17 ) 

where A, B and k are constants. Therefore, the field theory limit can be easily 
taken since we can let po infinity. 

In Table 1, we show the values of A, B and k for some values of the coupling 
constant go- The central charge becomes A + B at the held theory limit. The 
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calculated values of the central charge are shown as the function of the coupling 
constant go in hg.l. We also plot the central charge calculated for the sinc-Gordon 
theory by Itoyama and Moxhay [9]. As can be seen from the fig. 1, the two values 
of the central charge are quite different from each other. 

How can we interprete these differences ? The first possibility is that the two 
theories (sine-Gordon and massive Thirring models) are different from each other 
at the finite volume. We do not know whether this difference can show up as 
the central charge or not. However, a simple-minded physical intuition suggests 
that the central charge which should correspond to the heat capacity cannot be 
different if all the correlation functions of the two models are the same with each 
other. In this case, we should rather check the convergence of the perturbation 
expansions in Coleman’s proof of the equivalence between the sine-Gordon and 
the massive Thirring models since it crucially depends on the convergence of the 
expansions. For the negative values of the coupling constant, we do not know 
whether this convergence is already verified or not. 

The second possibility is that neither of the calculations are accurate enough 
to argue the difference between them. To this, we should comment on the ac¬ 
curacy of the present calculations. Since we have only the limited number of 
particles, we always face the criticism that the real nature (even though 1+1 
dimension) must be with the infinite number of particles. We have varied the 
number of particles from 1000 to 10000. It seems to us that the extracted central 
charge may well be reliable to within a few tens of percents. At least, we believe 
that the calculation must be rather reliable for the coupling constant around 
g 0 = — | where the extracted central charge is not very small. On the other 
hand, the present calculation may involve somewhat large errors for the coupling 
constant around or smaller than go = — | since the extracted central charge is 
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rather small. This is in contrast to the bound state problems [5,12-13] where 
there is some possibility of controlling the accuracy of the numerical calculations. 
However, the evaluation of the central charge involves rather complicated pro¬ 
cesses of extracting it since we have to obtain it from the term proportional to 
Y in the vacuum energy. Therefore, the error bars of the calculations we have 
shown in fig. 1 may well be still optimistic numbers. 

Concerning the central charge of the sine-Gordon model, we do not know 
whether the central charge predicted by Itoyama and Moxhay can be taken to 
be exact or not. Here, we only make a comment on the string hypothesis in the 
massive Thirring model when they employ the thermodynamic Bethe ansatz [14]. 
As discussed in ref. [5,11], the string picture in the massive Thirring model in the 
positive values of the coupling constant turns out to be invalid in the sense that 
they do not satisfy the PBC equations. However, in the negative values of the 
coupling constant, we do not know whether there is a string ^-like solution that 
satisfies the PBC equations. 

We would like to thank M. Hiramoto, C. Itoi, M. Kato and H. Mukaida for 
helpful discussions and comments. 
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Table 1 


90 

7r 

A 

B 

K 

c(go) 

-1/4 

0.941 

- 0.562 

25 

0.38 ± 0.09 

- 0.276 

1.06 

- 0.745 

16.5 

0.32 ± 0.05 

- 0.291 

1.06 

- 0.795 

16 

0.27 ± 0.06 

- 0.305 

0.901 

- 0.739 

25 

0.16 ± 0.07 

- 0.319 

0.854 

- 0.790 

30 

0.06 ± 0.04 

-1/3 

0.793 

- 0.879 

40 

- 0.09 ± 0.10 


We show the values of A, B and k for some coupling constant g 0 
together with the c(go) at the held theory limit. 
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